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Abstract — In  order  to  achieve  high  data  rate  digital  com¬ 
munications,  multiple-input/multiple-output  (MIMO)  techniques 
have  attracted  growing  interests  in  the  underwater  acoustic  com¬ 
munication  studies.  In  this  paper,  multichannel  combining  and 
decision  feedback  equalization  (MCC/DFE)  has  been  proposed 
for  underwater  acoustic  MIMO  channels.  In  order  to  overcome 
the  difficulties  introduced  by  the  fast  fluctuating  channel,  Doppler 
tracking  and  frequent  channel  estimation  are  performed.  Then 
time  reversal  combining  followed  by  a  single  channel  DFE  is  used 
to  demodulate  individual  symbol  sequences  transmitted  by  the 
multiple  element  source.  To  improve  the  performance,  successive 
interference  cancellation  is  also  incorporated  into  the  receiver 
structure. 

Using  data  from  the  Makai  experiment  conducted  around 
Kauai  Island,  HI,  2005,  we  have  shown  that  the  achievable  data 
rate  can  be  increased  up  to  4  times  using  the  same  bandwidth  as 
single  source  systems.  For  example,  32  kilobits/s  could  have  been 
achieved  by  simultaneous  transmission  of  four  4  kilosymbols/s 
4-phase  shift  keying  (QPSK)  symbol  sequences  when  both  the 
source  and  the  receiver  were  drifting  at  a  2  km  range  in  the 
ocean. 


I.  Introduction 

Underwater  acoustic  communications  can  provide  a  flexi¬ 
ble  and  cost-effective  way  to  exchange  information  between 
underwater  platforms.  However,  high  data  rate  underwater 
acoustic  communications  is  challenging  because  of  severe 
multi-path  spread,  time  varying  property  and  limited  band¬ 
width  of  the  channel  [1].  As  being  bandwidth  efficient  in  a 
rich  scattering  environment  [2]— [4],  multiple-input/multiple- 
output  (MIMO)  techniques  have  attracted  growing  interests 
in  the  underwater  acoustic  communication  studies.  Based  on 
the  multichannel  DFE  structure  proposed  by  Stojanovic  et  al 
[5]— [7],  various  receivers  for  MIMO  systems  were  developed 
with  a  focus  on  the  utilization  of  error  correcting  codes  [8]- 
[10].  In  [1 1]— [13],  the  time  reversal  based  MIMO  receivers 
were  developed  by  taking  advantage  of  spatial  diversity  of 
underwater  acoustic  channels. 

In  this  paper,  a  multichannel  combining  and  decision  feed¬ 
back  equalizer  (MCC/DFE)  has  been  proposed  for  underwater 
acoustic  MIMO  channels  based  on  the  receiver  for  single 
source  systems  presented  in  [14].  Using  data  obtained  from 
the  Makai  experiment  (MakaiEx)  [15],  [16]  conducted  around 
Kauai  Island,  HI,  2005,  we  have  shown  that  the  achievable  data 
rate  can  be  increased  2  to  4  times  using  the  same  bandwidth  as 


single  source  systems.  In  Section  II,  the  MIMO  system  model 
is  briefly  presented.  The  receiver  structure  is  presented  in 
Section  III.  The  experimental  results  are  shown  in  Section  IV. 
In  the  following  sections,  a  variable  with  "denotes  the  estimate 
of  the  variable,  c*  denotes  the  complex  conjugate  of  a  complex 
number  c.  a(n)*6(n)  denotes  the  convolution  of  two  sequences 
a(n)  and  b(n). 

II.  System  Model 

Consider  an  underwater  acoustic  MIMO  system  with  Nt 
transducers  and  Nr  hydrophones.  At  the  l- th  transducer  at  the 
source,  an  information  symbol  sequence  xi(n)  is  modulated 
to  carrier  frequency  fc  and  transmitted.  All  Nt  symbol 
sequences  from  Nt  transducers  are  independent  of  each  other 
but  with  the  same  symbol  rate  R  and  carrier  frequency  fc.  Let 
ym(t)  be  the  received  baseband  signal  at  the  m-th  hydrophone. 
The  effect  of  the  transmission  medium  between  the  l- th 
transducer  and  the  m-th  hydrophone  can  be  characterized 
by  a  time-varying  channel  impulse  response  (CIR)  function, 
Then  the  received  signal  on  the  m-th  hydrophone 
ym(t )  is  the  summation  of  all  symbol  sequences  distorted 
by  the  channel.  The  analog  waveform  ym(t)  is  sampled  at 
a  fractional  symbol  interval  to  provide  robustness  to  carrier 
phase  fluctuations  in  the  underwater  acoustic  channel  [5],  [17]. 
However,  for  notation  convenience,  symbol  spaced  signals  are 
used  throughout  the  paper.  Then, 

Nt 

Vmin)  =  ^V0i'm(n)  [xt (n)  *  hitm(n,fi)]  +vm(n),  (1) 

1=1 

where  ?/m(n)  is  the  discrete  time  representation  of  the  analog 
signal  ym(t ),  is  the  instantaneous  carrier  phase  offset 

associated  with  the  l- th  symbol  sequence,  and  vm(n)  repre¬ 
sents  the  ambient  noise.  /x),  0  <  y  <  L  —  1,  is  the 

discrete  time  baseband  CIR  function  where  L  is  the  duration  in 
symbols  for  all  the  channels.  /i/,m(n,  y)  includes  the  combined 
effects  of  transmitter/receiver  filters  and  the  CIR  function. 

The  challenge  in  underwater  acoustic  communications 
mainly  lies  in  the  highly  dispersive  and  time  varying  char¬ 
acteristics  of  the  channel.  At  high  data  rate  coherent  com¬ 
munications,  the  length  channel  L  is  usually  more  than  tens 
of  symbols.  Further,  multiple  hydrophones  are  employed  to 
achieve  desired  performance  for  even  single  source  underwater 
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Fig.  1.  The  MCC/DFE  structure  for  MIMO  systems. 


systems.  Therefore,  the  implementation  complexity  should  be 
taken  into  account  in  the  design  to  deal  with  multiple  highly 
dispersive  channels. 

III.  Receiver  Structure 

In  this  section,  a  low  complexity  MCC/DFE  structure  for 
underwater  acoustic  MIMO  systems  is  proposed  based  on  the 
receiver  in  [14].  The  MCC/DFE  combined  with  successive 
interference  cancellation  (SIC)  [17]  is  also  presented. 

A.  The  MCC/DFE  structure 

As  shown  in  Fig.  1,  the  receiver  consists  of  three  com¬ 
ponents:  Doppler  tracking  and  correction,  MIMO  channel 
estimation,  and  time  reversal  combining  and  single  channel 
DFE.  Doppler  tracking  and  correction  is  used  to  compensate 
for  Doppler  shift  and  any  other  linear  trend  in  the  carrier 
phase  offset.  MIMO  channel  estimation  is  performed  based  on 
the  Doppler  corrected  signals  and  multiple  estimated  symbol 
sequences.  To  demodulate  the  Z-th  symbol  sequence,  time 
reversal  combining  and  single  channel  DFE  are  followed. 

The  proposed  receiver  is  a  channel  estimation  based  struc¬ 
ture.  At  the  beginning  of  a  data  packet,  a  preamble,  or  a 
sequence  of  known  symbols,  is  used  to  perform  initial  channel 
and  Doppler  estimation  and  to  train  adaptively  the  DFE  tap 
weights.  After  the  preamble,  channel  and  Doppler  estimation 
are  frequently  updated.  The  most  recent  channel  estimate  is 


denoted  by  /x)  and  the  most  recent  Doppler  estimate 

is  denoted  by  /m(n).  The  major  parts  of  the  receiver  now  will 
be  discussed. 

1 )  Doppler  tracking  and  correction:  The  Doppler  estimate 
at  the  ra-th  hydrophone  is  obtained  by 


Na-1 


fm(n )  =  argmax 


^2  ym(n  -  p)(ym(n  -  p)ej2npfTs)* 

p= 0 


(2) 

where  ym(n)  =  Ya=i^i{u)  *  h,m(n,p)  and  Ts  =  1/R  is 

the  symbol  duration.  In  Eq.  2,  N/\  is  the  Doppler  observation 
block  size  in  symbols  and  <  f  <  +  - Sf , 

where  is  the  coarse  Doppler  estimate  and  5  f  is  the  Doppler 
search  range.  Since  Doppler  is  estimated  frequently,  is 
set  to  the  previous  Doppler  estimate.  The  Doppler  correction 
is  performed  by  offsetting  the  received  signal  2/m(n)  by  the 
estimated  Doppler  shift,  i.e.,  zm(n)  =  ^m(n)e_j2™^m(n)Ts, 
where  zm(n)  denotes  the  Doppler  corrected  signal. 


2)  MIMO  Channel  estimation:  Assuming  carrier  phase 
offset  6^m(n)  can  be  removed  completely  by  the  Doppler 
correction,  we  have 


nt 

zm ( n )  ="22x1  (n)  *  hi,m [n,  p)  +  v'm (n) ,  (3) 

1=1 


Fig.  2.  The  MCC/DFE  structure  with  SIC. 


where  vfm(n)  is  the  noise  term  after  Doppler  correction. 
The  NtNr  channel  estimates  h^m(n, /a)  can  be  obtained 
jointly  based  on  the  Doppler  corrected  signal  zm(n)  and  the 
previously  detected  symbols  xi(n)  or  the  known  symbols 
xi  (n).  Various  least  squares  algorithms  can  be  used  for  channel 
estimation.  In  this  paper,  the  iterative  least  squares  QR  (LSQR) 
algorithm  is  used  [18].  The  channel  estimation  block  size  is 
chosen  to  be  three  times  the  total  channel  taps  to  be  estimated, 
i.e.,  TVo  =  3NtL. 

3)  Time  reversal  combining  and  single  channel  DFE  for 
the  l-th  symbol  sequence:  To  demodulate  the  Z-th  symbol 
sequence,  the  other  symbol  sequences  can  be  treated  as  noise. 
The  time  reversal  combining  and  the  single  DFE  can  be  used 
to  demodulate  individual  symbol  sequence. 

Time  reversal  combining  uses  —/f))*  to  match- 

filter  the  Doppler-corrected  signals  on  each  channel  zm(n) 
and  then  combines  the  results  [19]— [21].  The  output  of  time 
reversal  combining  is 

Nr 

n{n)  =  Y  (hm(n,  -m))*  *  Zm(n) 

m= 1  '  ' 

=  xt{n )  *qi(n,n )  +wi(n), 
where  wi(n)  is  the  noise  component, 

Nr 

wi{n)  =  Y  -fi))*  *  (v'm(n)),  (5) 

m=  1 

and  qi(n,l)  is  the  effective  CIR  function,  or  the  q-function 
[22],  between  the  l-th  transducer  and  the  receiver: 

Nr 

qi(n,l)  =  ^  /'))*  *  hi,m(n,(i).  (6) 

m=  1 


A  single  channel  DFE  with  joint  phase  tracking  [6]  is  used 
to  equalize  the  residual  inter-symbol  interference  in  77(71).  The 
exponentially  weighted  recursive  least-squares  (RLS)  algo¬ 
rithm  is  used  to  update  the  equalizer  tap  weights.  The  residual 
carrier  phase  offset  in  77(71)  is  compensated  for  by  a  second 
order  phase  locked  loop  (PLL)  embedded  in  the  adaptive 
channel  equalizer.  The  phase  correction  based  on  the  PLL 
output  is  implemented  at  the  input  to  the  DFE  feedforward 
filter.  The  soft  output  signal-to-noise  ratio  (SNR),  p ,  of  the 
single  channel  DFE  is  used  as  a  performance  metric  in  this 
paper. 

B.  The  MCC/DFE  structure  with  SIC 

In  MIMO  systems,  multiple  symbol  sequences  are  simul¬ 
taneously  transmitted  and  each  symbol  sequence  causes  in¬ 
terference  in  the  demodulation  of  other  data  streams.  In  other 
words,  co-channel  interference  exists  because  multiple  symbol 
sequences  share  the  use  of  the  same  channel  and  frequency 
band.  In  order  to  mitigate  the  co-channel  interference,  SIC 
[17]  can  be  incorporated  into  the  MCC/DFE  structure.  The 
symbol  sequences  are  demodulated  in  the  order  of  the  soft 
output  SNR  of  the  single  channel  DFE.  As  shown  in  Fig.  2, 
the  strongest  symbol  sequence,  the  l\-th  symbol  sequence,  is 
to  be  demodulated  first.  After  the  Z^-th  symbol  sequence  is 
demodulated,  interference  resulted  from  it  will  be  removed. 
That  is, 

zm(n)  =  4r1}(n)  -  %u(n)  *  hh,m(n,fl),  (7) 

where  Zm{n)  denotes  the  Doppler  corrected  signal  with  the 
interference  removed  from  the  strongest  i  symbol  sequences. 
Then  the  /i+i-th  symbol  sequence  is  demodulated  based  on 

Zm(n). 


C.  Comparison  with  existing  acoustic  MIMO  receivers 

In  the  literature,  existing  acoustic  MIMO  receivers  include: 

(1)  multichannel  DFE  based  MIMO  receivers  in  [7]— [10]  and 

(2)  the  time  reversal  MIMO  receivers  in  [1 1] — [13].  In  the 
multichannel  DFE  based  MIMO  receivers,  feedforward  filters 
are  applied  to  the  individual  channels  and  their  outputs  are 
combined  prior  to  the  feedback  filter  as  in  the  multichannel 
DFE  [5].  Phase  synchronization  at  the  individual  channels  is 
optimized  jointly  with  the  equalizer  tap  weights.  The  number 
of  adaptive  feedforward  taps  increases  with  the  number  of 
hydrophones.  Unlike  the  multichannel  DFE  based  MIMO 
receivers,  the  proposed  receiver  uses  a  single  channel  DFE 
after  time  reversal  combining  for  each  symbol  sequence. 

An  advantage  of  the  proposed  receiver  structure  is  its  low 
complexity.  For  the  demodulation  of  each  symbol  sequence, 
the  complexity  of  a  multichannel  DFE  based  MIMO  receiver 
increases  at  least  as  the  square  of  the  number  of  hydrophones 
if  RLS  algorithms  are  used  for  a  fast  tracking  capability  [17]. 
Since  time  reversal  combining  collapses  multiple  channels  into 
a  single  channel,  the  complexity  of  the  successive  DFE  re¬ 
mains  unchanged  when  the  number  of  hydrophones  increases. 

Although  also  time  reversal  based,  the  proposed  receiver  has 
a  different  structure  than  the  time  reversal  MIMO  receiver  in 
[1 1]— [13]  where  multichannel  combining  is  performed  based 
on  channel  probes  or  the  known  symbols  at  the  beginning  of 
the  data  packet.  Phase  tracking  or  Doppler  tracking  usually 
is  performed  after  time  reversal  combining.  Compared  with 
the  referenced  time  reversal  MIMO  receivers,  the  proposed 
receiver  performs  continuous  Doppler  tracking  and  channel 
estimation  to  overcome  fast  fluctuations  which  occur  over  the 
duration  of  a  data  packet. 

IV.  Experimental  Results 

MakaiEx  was  conducted  from  Sept.  15  to  Oct.  2,  2005,  west 
of  Kauai,  HI,  to  study  high-frequency  underwater  acoustic 
communicatons  [15].  During  MakaiEx,  the  10  element  verc- 
tical  MIMO  source  provided  by  Space  and  Naval  Warfare 
Systems  Center  and  the  8  element  ACDS  receiving  array 
provided  by  Naval  Research  Laboratory  (NRL)  were  deployed 
three  separate  times  [16].  In  this  paper,  the  data  from  the 
second  MIMO  deployment  on  Sept  26,  2005  will  be  analyzed. 
The  8  element  ACDS  receiving  array  was  deployed  by  the  R/V 
Kilo  Moana  and  set  in  a  free  drift  mode.  The  top  element  of 
the  ACDS  array  was  about  20  m  below  the  sea  surface  and  the 
spacing  was  about  2  m.  The  sampling  frequency  of  the  ACDS 
array  was  160  kHz.  The  10  element  MIMO  source  was  hang 
from  the  deck  of  the  R/V  Kilo  Moana.  The  spacing  of  the 
source  element  was  2  m  with  the  top  element  about  20  m 
below  the  sea  surface.  The  source  power  level  of  each  source 
element  is  190  dB  re  1  pPa  at  1  m.  Once  deployed,  the  R/V 
Kilo  Moana  maintained  roughly  a  2  km  separation  with  the 
ACDS  array.  The  water  depth  of  the  experimental  site  was 
about  100  m. 

The  carrier  frequency  of  the  analyzed  communication  signal 
is  fc  =  37.5  kHz  and  the  symbol  rate  is  R  =  4  kilosymbols/s. 
The  square-root  raised  cosine  shaping  filter  is  used  with 
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Fig.  3.  The  formation  of  the  data  packet. 


an  excess  bandwidth  [17]  of  75%.  The  communication  data 
were  in  a  form  of  packets.  A  1248  symbol  long  preamble 
preceded  the  data  packet.  Then  each  800  symbol  data  block 
was  protected  by  a  224  symbol  prefix  and  a  224  symbol 
appendix.  The  prefix  and  appendix  were  used  to  periodically 
retrain  the  equalizer.  After  the  preamble,  7  data  blocks  were 
transmitted  continuously.  Note  the  7-th  block  had  a  variable 
length  less  than  800  symbols.  The  total  length  of  the  packet 
was  about  2.5  s.  The  data  packet  is  illustrated  in  Fig.  3. 
Three  types  of  source  options  were  used,  i.e.,  1  transducer, 
2  transducers,  and  4  transducers,  to  transmit  binary  phase 
shift  keying  (BPSK)  and  4  phase  shift  keying  (QPSK)  signals. 
A  total  of  6  packets  are  shown  to  demonstrate  the  receiver 
performance. 
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Fig.  4.  The  CIR  function  on  the  top  receiving  element  of  the  ACDS  array.  The 
depths  of  the  source  elements  are  22  m,  26  m,  32  m,  and  38  m,  respectively. 
Color  scale  has  a  20  dB  dynamic  range. 


As  mentioned,  fractional  spaced  sampling  is  used  in  the 
receiver  and  the  oversampling  rate  is  K  =  4.  The  estimated 
length  of  the  CIR  function  is  25  ms,  or  L  =  100  symbols, 
for  the  single  element  source  packets.  It  is  set  as  L  =  50 
symbols  for  the  multiple  element  source  packets.  The  channel 
estimation  block  size  No  and  Doppler  observation  block  size 
A/a  are  both  set  as  3NtL.  The  channel  estimation  update 
interval  is  chosen  as  N  =  100  symbols.  At  the  beginning  of 
the  packet,  1248  symbols  are  used  to  conduct  initial  channel 
estimation,  Doppler  tracking,  and  DFE  tap  weight  training. 
The  Doppler  search  range  is  Sf  =  1.6  Hz.  The  number  of  the 
feedforward  taps  is  KNff  =  40  symbols  for  the  fractionally 
spaced  DFE  [6]  where  Nff  =  10  is  the  feedforward  filter 
span  in  symbols.  The  number  of  the  feedback  taps  is  Nfb  =  2 
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TABLE  I 

Receiver  parameters. 


Parameters 

Description 

Value 

fs 

Sampling  rate 

160  kHz 

fc 

Carrier  frequency 

37.5  kHz 

R 

Symbol  rate 

4  kHz 

fEB 

Excess  bandwidth  of  the  square-root  raised  cosine  filter 

3  kHz 

K 

Oversampling  factor 

4 

Nr 

Total  number  of  the  hydrophone  channels 

8 

Npreamble 

Size  of  the  preamble 

1000  symbols 

L 

Length  of  the  CIR  function 

50  or  100  symbols 

N0 

Channel  estimation  block  size 

2 JNtN  symbols 

N 

Channel  estimation  update  interval 

100  symbols 

na 

Doppler  observation  block  size 

3 Nt-L  symbols 

V 

Doppler  search  range 

1.6  Hz 

Nff 

Feedforward  filter  span  in  symbols 

10  symbols 

Nfb 

Feedback  filter  tap  number 

2  symbols 

Kh 

Proportional  tracking  constant  in  PLL 

0.0002  I 

Kf2 

Integral  tracking  constant  in  PLL  ~  ' 

0.0002  ‘ 

X 

RLS  forgetting  factor  in  the  DLE 

because  the  feedback  filter  is  applied  to  a  symbol  spaced 
sequence.  The  RLS  forgetting  factor  A  in  the  DFE  is  chosen 
as  0.999.  In  the  PLL  embedded  in  the  DFE,  the  proportional 
tracking  constant  Kfx  and  the  integral  tracking  constant  Kf2 
are  both  set  as  0.0002.  The  main  parameters  are  listed  in 
Table  I. 

Fig.  4  shows  the  selected  CIR  functions  for  the  4  transducer 
BPSK  packet.  As  shown,  the  CIR  functions  exhibits  fast  fluc¬ 
tuation  during  the  2.5  s  packet.  Such  fast  fluctuating  channels 
require  frequent  channel  estimation.  The  other  detrimental 
effect  of  the  underwater  acoustic  channel  is  the  fast  phase 
fluctuation.  In  the  MCC/DFE  structure,  it  is  treated  by  the 
Doppler  tracking  and  correction.  Fig.  5(a)  shows  the  Doppler 
estimate  at  the  four  receiving  elements  for  the  4  transducer 
BPSK  packet.  As  shown,  the  Doppler  estimates  of  different 
elements  are  largely  correlated.  The  instantaneous  Doppler 
estimates  can  be  more  than  4  Hz  for  the  drifting  source  and 
receiver.  When  the  Doppler  correction  has  been  made,  the 
residual  phase  estimate  is  slowly  variant  as  shown  in  Fig.  5(b). 

The  receiver  performance  is  listed  in  Table  II.  For  all 
the  packets,  the  two  receivers  are  able  to  track  the  channel. 
The  maximum  data  rate  is  32  kilobits/s  achieved  by  the  4 
transducer  QPSK  packet.  The  overall  bit  error  rate  (BER)  of 
the  packet  is  0.06  for  the  receiver  with  SIC.  SIC  can  improve 
the  output  SNR  up  to  3  dB.  However,  for  the  4  transducer 
QPSK  packet,  the  improvement  is  minimum  (less  than  1  dB 
for  all  symbol  sequences).  The  demodulation  order  in  SIC 
is  preset  according  to  the  soft  output  SNR  in  the  MCC/DEF 
receiver. 

V.  Conclusion 

In  this  paper,  a  MCC/DFE  structure  has  been  proposed  for 
underwater  acoustic  MIMO  channels.  In  order  to  overcome  the 
difficulties  introduced  by  the  fast  fluctuating  channel,  Doppler 
tracking  and  frequent  channel  estimation  are  performed.  Then 
time  reversal  combining  followed  by  a  single  channel  DFE  is 


used  to  demodulate  individual  symbol  sequences  transmitted 
by  the  multiple  element  source.  SIC  is  also  incorporated 
into  the  structure  to  improve  the  receiver  performance.  Using 
experimental  data  from  MakaiEx,  we  have  shown  that  the 
achievable  data  rate  can  be  increased  2  to  4  times  using  the 
same  bandwidth  as  single  source  systems. 
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